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Abstract 

As an ensemble scheme of solid-state NMR quantum computers the extension of Kane's many- 
qubits silicon scheme based on the array of 31 P donor atoms are spaced lengthwise of the strip gates 
is considered. The possible planar topology of such ensemble quantum computer is suggested. The 
estimation of the output NMR signal was performed and it was shown that for the number N > 10 5 
of ensemble elements involving L ~ 10 3 qubits each, the standard NMR methods are usable. 

As main mechanisms of decoherence for low temperature (< O.IK), the adiabatic processes 
of random modulation of qubit resonance frequency determined by secular part of nuclear spin 
hyperfine interaction with electron magnetic moment of basic atom and dipole-dipole interaction 
with nuclear moments of neighboring impurity atoms was considered, It was made estimations 
of allowed concentrations of magnetic impurities and of spin temperature whereby the required 
decoherence suppression is obtained. Semiclassical decoherence model of two qubit entangled states 
is also presented. 

As another variant of the solid-state ensemble quantum computer, the gateless architecture of 
cellular-automaton with antiferromagnetically ordered electron spins is also discussed here. 
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Atomic nuclei with spin quantum number 1 = 1/2 are the natural candidates for qubits in quantum 
r S^ computers. The early approach to NMR quantum computers was suggested in 1997 E] and then 
confirmed in experiments jHJ Q] . In this approach several diamagnetic organic liquids whose individual 
molecules, having a number of interacted non-equivalent nuclear spins-qubits with 1 = 1/2 and being 
nearly independent on one another where used. They act in parallel as an ensemble of almost indepen- 
dent quantum molecules-microcomputers. In so doing the nuclear spins of an individual molecule are 
described by mixed state density matrix of reduced quantum ensemble. Initialization of the nuclear 
spin states in this case means the transformation of mixed state into so called, effective or pseudo-pure 
state UJHHI3]. 

The access to individual qubits in a liquid sample is replaced by simultaneous access to related 
qubits in all molecules of a bulk ensemble. Computers of this type are called bulk-ensemble quantum 
computers. The liquid-based quantum computer can operate at room temperature. For control and 
measurements of qubit states the standard NMR technique is used. 

The principle one-coil scheme of experiment is shown in Fig. 1. The sample is placed in the 
constant external magnetic field B and in the alternating (say, linearly polarized) field b(t), produced 
by RF voltage K,(t): 

B(t) = B+b(t) = Bk+2bcos(ujt + tp)\, (1) 
where i and k are unit vectors along the axes x and z. 
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Fig. 1. The principle one-coil scheme of NMR measurement. 

Let the sample represent an ensemble of N molecules-microcomputers with L qubits each at 
temperature T = 300 K, in the external magnetic fields B = 1 — 10 T. The resonance nuclear spin 
frequency is uj a /2it ~ jiB/2n < 150 MHz, 71 is gyromagnetic ratio of nuclear spin (71 ~ 7n = 
95.8radMHz/T), hto A /kT < 10~ 5 . 

The output oscillating voltage V(t) is 

V(t) = QKd${t)/dt = noQKAdM x (t)/dt, (2) 
where = j A [i$M x (t)dydz is magnetic flux produced by resonant spins in the coil (fio = Air ■ 

1Q -1 T 2 

cm 3 /J), L s = hq(KA) 2 /V s is solenoid inductance of the resonance counter, V s is volume of 
the solenoid, K is the number and A is area of coil turns, Q = R/{uaL s ) > 10 2 is the quality 
factor of resonance counter for parallel connected resistance R (Fig. 1). For resonance condition 

u = uj A = (L s cr 1 / 2 . 

The maximum nuclear spin read-out magnetization M xmax (the liquid sample is considered here 
to be a continuous medium and to have volume V ~ V s ) at optimum resonance condition is defined 
by the amplitude of RF field b = 1/(71 y^TnTj|i) [H] (see also (20) below): 

M xmax = M zmV /T ±I /T|| I /2 w ■ (N/V s ) ■ e(L)/2, (3) 

where M zm is maximum equilibrium nuclear magnetization, T±j and Tp are effective transverse and 
longitudinal relaxation times, N is number of resonant nuclear spins (one in a molecule) in volume 
V s . Parameter e(L) is the maximum probability of the full nuclear polarization in pseudo-pure state 
Pi = 1 It may be estimated by the difference of equilibrium population between the lowest and 
the highest energy states. For nearly homonuclear L spin system j7] it is: 



e(L) 



ex.p(LhtJA/2kT) — exp(—LhujA/2kT) 2sm\\{LTiuj A /2kT) 



(exp(hu A /2kT) + exp(-7kj A /2£;T)) L 2 L cosh L (^ A /2A;T) 



(4) 



In the high temperature limit huo\/(kT) C 1 we have e(L) = L2~ L %u A / '(kT), that is, the signal 
amplitude exponentially drops with the number of qubits, but it does not drop for huj A / (kT) S> 1 
when e(L) = 1 (the pure ground nuclear spin quantum state). 

The maximum NMR signal intensity S is defined by amplitude 



S = l^maxl = (w/4)QKA(N/V B )7inu A e(L), 
where the product KA can also be expressed as 

KA = (L s F s //x ) 1/2 = (PV s /(W)^a)) 1/2 . 
For the root-mean square noise voltage in the measurement circuit we write 

V N = VAkTRAu, 



(5) 



(6) 



(7) 
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where the amplifier bandwidth is as a rule Av ~ 1 Hz. 
So for signal to noise ratio we obtain 

(S/N, = |K TO |/ yN ,l^|g ^ W ~ 

~ 0.2yJ(Q/V s ) ■ {hLO A /kT)Ne(L) • 10~ 9 , (here F s in cm 3 ). (8) 

For example, for two qubits molecules (L = 2), using, e(L) = hu>\/(2kT) ~ 10 -5 , we can make an 
estimation 

(S/N) ~ (Q/V S ) 1/2 N • 10~ 16 . (9) 

Thus, to keep the value (S/N) > 1, the number of resonant nuclear spins for two qubit liquid ensemble 
at room temperature, V s ~ 1 cm 3 and Q ~ 10 3 is bound to be N > 10 16 . 

In the case of paramagnetic liquids one would expect that the number of polarized nuclei may 
be increased with dynamic polarization (say, Overhauser effect). Assuming electron and nuclear 
gyromagnetic ratio Je/"fl ~ 10 3 we obtain that in the probability e(L) for a L— qubits single state the 
value hu>A/(kT) in (JHJ) should be replaced by W 3 hiOA/{kT). Therefore, for the same value e(L) and 
number of molecules N, the allowed number of qubits L approximately will be estimated from 

L2~ L > 1(T 3 , (10) 

whence it follows that L < 12 qubits. 

An additional increase of read-out NMR signal may be obtained in paramagnetic liquids using the 
ENDOR technique. It is generally believed that for the liquid bulk-ensemble quantum computers a 
limiting value is L < 20 — 30 0. 

There are five basic criteria for realization of a large-scale NMR quantum computer, which can 
outperform all traditional classical computers jHJ: 

1. For any physical system, which presents large-scale quantum register, the necessary number of 
qubits in quantum register must be L > 10 3 . 

One such example of this register is solid-state homonuclear system, in which nuclear spin 
containing identical atoms are housed at regular intervals in a natural or an artificial solid-state 
structure. 

2. There is a need to provide the conditions for preparation of initial basic quantum register state. 
For a many-qubit solid-state NMR quantum computer the quantum register state initializing 
can be obtained by going to extra-low nuclear spin temperature (< lmK at fields of order of 
several tesla). 

3. The decoherence time of qubit states should be at least up to 10 4 times longer than the 'clock 
time', that is value of order of several seconds for NMR quantum computers. The decoherence 
suppression is one of the important problems in realization of a large-scale quantum computers. 

4. There is a need to perform during a decoherence time a set of quantum logic operations deter- 
mined by a logic unitary transformation. This set should contain certain set of the one-qubit 
and two-qubits operations are shielded from random errors. The electromagnetic pulses that 
control the quantum operation should be performed with an accuracy of better than 10 -4 -10 -5 . 

5. There is a need to provide accurate and sensitive read-out measurements of the qubit states. This 
is another of the important and hard problems. 

The design of solid-state NMR quantum computers was proposed by B. Kane in [HI QH]. It was 
suggested to use a semiconductor MOS structure on a 28 Si spinless substrate, in a near-surface layer 
whose stable phosphorus isotopes 31 P, acting as donors, are implanted in the form of a regular chain. 
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These donors have a nuclear spin 1 = 1/2 and substitute for silicon atoms at the lattice sites, producing 
shallow impurity states. The number of donors or the qubit number L in such a quasi-one-dimensional 
artificial 'molecule ' may be arbitrary large. It is suggested an individual nuclear spin-qubits electrical 
control and measurement of qubit states through the use of special gate structures. The experimental 
implementation of Kane's scheme is undertaken now in Australian Centre for Quantum Computer 
Technology [TT1 IP2]. 

However, there are four essential difficulties in implementing this quantum computer: 

1. First of all, signal from the spin of an individual atom is very small and high sensitive single-spin 
measurements are required. 

2. For initialization of nuclear spin states it is required to use very low nuclear spin temperature 



3. It is required to use regular donors and gates arrangement with high precision in nanometer 



4. It is necessary to suppress the decoherence of quantum states defined by fluctuations of gate 



As an alternative, we proposed the variant of an ensemble silicon-based quantum computer |13l I14j . 
One would expect that with the ensemble approach, where many independent 'molecules' of Kane's 
type work simultaneously, the measurements would be greatly simplified. Here we will give some 
further development of this scheme. 

1 The silicon structure with regular system of strip gates 

In this case, unlike the structure suggested in 9 , gates A and J form a chain of narrow (I a ~ 10 nm) 
and long strips along which donor atoms at l y distant from each other are placed (Fig. 2). Thus, they 
form a regular structure of the planar silicon topology type. 
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Fig. 2. The structure of two qubit cells for three ensemble component. 
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The separation between neighboring donor atoms in Si, as in Kane's scheme, must be l x < 20 nm. 
In this case the interqubit interaction is controlled by gates J. The depth of donor d is ~ 20 nm. For 
l y 3> l x the exchange spin interaction between electrons of donor atoms disposed along the strip gates 
(y— axis) is negligibly small. Hence, such a system breaks down into an ensemble of near-independent 
Kane's artificial 'molecule', whose electronic spins at temperature T < 0.1K are initially fully aligned 
with the field of several Tesla (^ e TiB/kT ^ 1). As in case of liquids, the nuclear spin states of 
individual Kane's chain- 'molecule' will be described by density matrix of reduced quantum ensemble. 
Access to individual qubits will be replaced by simultaneous access to related qubits in all 'molecules' 
of ensemble. 

The linear qubit density in the artificial 'molecules' is ~ 50 qubits on micrometer. For the real- 
ization of considered structure, as well as of the Kane's scheme, the nanotechnology with resolution of 
the order of ~ 1 nm is also needed. 

For the initializing of all nuclear spin-qubit quantum states (fully polarized nuclear spins) there is 
a need to attain, for the time being, nuclear spin temperature T < 10 _3 K. An output signal in this 
system, as in liquids, will be proportional to the number of 'molecules' or donor atoms N (component 
number of our ensemble) in the chain along axis y. In the following the lower value of N will be 
estimated. 



2 The states of insulated donor atoms in magnetic fields 

The electron-nuclear spin Hamiltonian for a donor atom 31 P has the form 

H = 7 c nBS-7i^BI+AIS, (11) 

four energy levels of which are given by the well-known Breit-Rabi formula. For / = 1/2, 5 = 1/2 
(the z— axis is parallel to B) this formula is written as 

A 



E(F,m F ) = ----/ihBm F - (-l) F sign(l + m F X) — y/l + 2m F A + X 2 , (12) 

where constant of hyperfine interaction A/(2irh) = 116 MHz 15 , X = (j c ± ji)hB /A ~ ~f c TiB/A 3> 1, 
F = I± 1/2 = 1,0, and m F = M + m = ±1,0, if F = 1 or m F = 0, if F = (Here M = ±1/2 and 
m = ±1/2 are z— projections of electron and nuclear spins accordingly). The energy level scheme is 
shown in Fig. 3. For the energy of the ground spin state, F = and m F = 0, hence, we obtain 



E(0, 0) = -A/A - (A/2)Vl + X 2 . (13) 

For the next, excited energy state, F = 1, m F = —1 we have 

E{l,-l)=A/A-( 7e - 11 )hB/2. (14) 

Thus, the energy difference between the two lower states of the nuclear spin (the resonant qubit 
frequency), that interacts with an electron, whose state remains unchanged, is described in simple 
terms (7 e > 71 for X w j c hB/A > 1): 

hu£ = E(l, -1) - E(0, 0) = A/2 + (71 - 7 e)h5/2 + jy/l + X 2 « 

A A 2 
-fl%B + - 



2 AjehB' 

A 

2 A le hB' 



hcjl = £(1,1) -£(1,0) ~ -71^ + 7T + 77^- (15) 



For 31 P donor atoms 7 e /7i = 1.62 • 10 3 , j c = 176.08 radGHz/T, ji = 1.13tn = 108radMHz/T. In 
magnetic field B = IT: ^/2vr = 75 MHz, uj^/2tt = 41 MHz. 
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Fig. 3. Energy levels of an individual donor atom in magnetic field. 



The frequencies los, wb, ^c, W D are in microwave, - in the RF ranges of frequencies. The 
transitions with frequencies lus in the first approximation are forbidden. 
The states |jP, mp) in M,m basis are 



11,-1) 
|i,o) 

|0,0) 
a 



|l/2,l/2> , 
1-1/2,-1/2), 

(1 - a) 1 ' 2 |l/2, -1/2) + a 1 ' 2 |-l/2, 1/2) , 
(l-a)V2|_ 1/2 , 1/2)- ^11/2, -1/2), 



1 



1 



1/(4X 2 ) < 1. 



(16) 



The diagonal matrix elements of nuclear magnetization M z per one donor atom for two lower energy 
states will be determined by 



(0,0|M Z |0,0) 
(1,-1|M Z |1,-1) 



(0,0| I z |0,0)7ift 

(l,-l|iz|l 



X 



VT+x 2 ' 
i) iih = -nh/2 



(17) 



The probabilities of the L— qubit lowest and highest energy fully filling states for the same electron 
spin state M = —1/2, which correspond, as noted above, to the maximum probability of the nuclear 
polarization in pseudo-pure state, are: 



/(1,-1) 



p L (0,0) 



exjp(—LTitu~^/2kT) 



(exp(hujX/2kT) + exp(-7k4/2/cT)) L ' 

exp(LhujX/2kT) 
(exp(hu;£/2kT) + exp(-huj£/2kT)) L ' 



(18) 



The possible maximum nuclear magnetization M zm (the populations of states |1, 1) and 1 1.0) is negli- 
gible for ws,wb,^c > w a) i s 



M 7 - 



7ih/2 ■ (N/V c ) 



X 



ex.p(Lhoj~X /2kT) 



y/l + X 2 (exp(hu£/2kT) + exp(-^/2fcT)) i 
exp(-L^/2fcT) \- ll K, 2 .(N/V c )e(L). 



(exp(huj£/2kT) +exp(-^X/2A:T)) L 

For Lhuj£/2kT <C 1 and 1 > 1 we obtain (compare with (pj)) 

M zm « 7 in/2 • (JV/K) • 2~ L L(hiuX/kT). 

But for very low temperatures {%uj~£/2kT 3> 1) we have the /u// nuclear polarization M z 
(N/V c ) and e(L) = 1. 



(19) 

(20) 
7tty2 . 
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3 The gain effect for NMR signal 



Transitions between two lower states are induced by a RF magnetic field, applied at a frequency 
resonant ut- The Rabi resonance frequency Q, which is defined by matrix elements of spin interaction 
Hamiltonian with the external RF field b(t) 

flrf (*) = (lcS x - Ti^)^x(t), b x (t) = 26cos(u;+i) (21) 

can be found from 

n = 11 b eS (X) =2 |<0, 01^(0)11,-1)1 /ft. (22) 
For the amplitude of effective RF field, acting on nuclear spin, b e s(X) we obtain 

b cS (X) = b (a^dehi) + (1 - a) 1 / 2 ) , (23) 

where b is the amplitude of circularly polarized field component. 

The Rabi frequency has the maximum value for X = (a = 1/2) and monotonically reduces to 
value for the insulated nuclear spin (a 0), ^fib e s(X 3> 1) = 716. From the rate of quantum operation 
standpoint it is desirable to operate in relatively weak fields jl()j . at which r y e /ji 3> X ~ j e hB/A 3> 1 
or 3.5 T > B > 3.9 • lO" 3 T. 

In this case from (|23j) we will obtain 

b cS = (1 + rj)b » 6, (24) 

where r\ = A/(2~fih~B) 3> 1 is the gain factor. Under these conditions RF field operates through 
the transverse component of electronic polarization. For magnetic fields B = IT we have the value 
b e g = 4.4 • b, and for B = 0.01 T we have the value b c s = 338 • b. The gain effect involves an increase 
of NMR signal and Rabi frequency. This effect was indicated previously by K. Valiev in [lfij . 

In the pulse technique this effect makes it possible to decrease the length of pulse and along with it 
the times of logic operation performing. Moreover, the computer operations, owing to this effect, can 
be performed at lower RF fields. At last, it permits to reduce the RF field influence on the operation 
of neighboring semiconductor devices. 

To describe the nuclear dynamics for the two low-lying level systems being discussed (X ^> 1), we 
can write the following Bloch-type equation with only two effective relaxation times: 

^ = 7l[ MxB „]-^i- < M -- M -> k (25) 

where i, j, k are orthogonal unit vectors (Fig. 1), M zm is defined as (|T§|). 

Bcflf = (wA/7i) k +2freff cos(u;£)i. (26) 

It follows from it that the value of maximum nuclear read-out magnetization in NMR signal has 
here for 6 e ff(X) = 1/ (71 ^JTjjT^i ) , that is again for M xmax = M zm y / Tj_i/Tj| I /2. Hence, the read-out 
NMR signal can not be increased through the gain effect over its maximum value, that corresponds 
to M zm ,/r ±I /T|| I /2. 



4 The signal to noise ratio for an ensemble silicon quantum com- 
puter 

For the realization of an ensemble silicon quantum register we propose a variant of planar scheme |17j . 
that, as an example, contains n-p in parallel acting identical blocks, each has A^o in parallel connected 
L— qubit Kane's linear 'molecules'. This scheme is schematically depicted in Fig. 4. 

Let the sample be the silicon ( 28 Si) plate of thickness 0.1 cm. For the full number of computers- 
' molecules' in ensemble N = p- Nq ■ n, the volume of sample and also of solenoid is V s ~ S ■ l x ■ l y ■ L ■ N 
(the filling factor is assumed for simplicity to be one). 
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Fig. 4. The scheme of the proposed planar silicon topology with p ■ n in parallel connected blocks of 
the ensemble L— qubit quantum computers (the connections are not shown here). The broad and 

narrow lines denote the A and J gates. 



The read-out signal from such ensemble in parallel acting chains, as distinct from liquid prototype, 
for full nuclear polarization or, what is the same, for nuclear spin temperatures T\ < 10 -3 K has instead 
of the small factor in intensity of the NMR signal of type s(L) = 2~ L L ■ TiujxI (kT) the factor e(L) = 1. 
The NMR signal from our sample within a non-essential factor is the same as from macroscopic sample 
(see Appendix A.l). Therefore, with the expressions (@J and (jHJ), %u&/(kTi) < 1 (Tj < lmK) and 
e = 1 we will obtain as an estimation for maximum signal to noise ratio 

(S/N) « ^Qhu A /(kTV s ) ■ N ■ 1(T 9 « ^QN/iSlJyL) • 1(T 10 . (27) 

It is believed that for low temperatures Q ~ 10 6 . The effective volume of one 'molecule' for l x = 20 nm, 
l y = 50 nm, L = 10 3 , V s = 5l^l y L = 10" 9 cm 3 we receive that the read-out signal in our scheme may be 
available for standard NMR technique, if the number of 'molecules' in ensemble is of about iV > 10 5 . 
So high-sensitive devices for measurement of individual spin-states are not needed. 

To estimate the values n, p let us consider the square plate with 50Nqp = 20 • 10 3 n and iVo = 100. 
As a result, we receive n ~ 16 and p ~ 63. The area of the structure without passive regions is ~ 315 
x315 fim 2 . This size is sufficiently small for sample to be housed in the split between the magnet poles 
of a standard NMR spectrometer. Real plate may have considerably more area and correspondingly 
more number of 'molecules' N. 

For implementation of two-qubit logic operation it is required the controlled by gates J interqubit 
indirect interaction with characteristic frequency vj ~ 100 kHz <C wa/2vt ~ 100 MHz. To bring about 
fault-tolerant quantum computations on large-scale quantum computers the relative error for single 
logic operation must not be more than ~ 10 _5 |S]. Hence it follows that a resolution bound of the 
NMR spectrometer must be of the order of ~ 100 kHz • 10~ 5 ~ 1Hz, that is consistent with the usual 
requirements. It is significant that such high precision is needed only for performing the logic quantum 
operation, but it is not needed for read-out measurements. 

The read-out signal may be more increased by means of an electron-nuclear double resonance 
(ENDOR) methods JH] of observing the electron resonance at transition with frequencies ui-q and uq 
(Fig. 3). 

Consequently, by the use of standard NMR and additional of ENDOR techniques the first main 
difficulty of Kane's scheme can be overcome. 
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5 The cooling of nuclear spin system and nuclear state initialization 
by means of dynamic polarization 

The electron and nuclear longitudinal relaxation times for the allowed transitions in four energy level 
system of phosphorus doped silicon have been extensively investigated experimentally in |18| 119] . 
For the allowed transitions with frequency wb and loq (Fig. 3) electron longitudinal relaxation times 
r ||B ~ r ||c a t l° w temperatures were found to be exceedingly long. They are of the order of one 
hour at T = 1.25 K, B ~ 0.3 T, are independent of phosphorus concentration below C ~ 10 16 cm -3 
(mean distance between phosphorus atoms is of the order of 45 nm) and are approximately inversely 
proportional to the lattice temperature T. The nuclear longitudinal relaxation time Xli (the frequency 
c</£) were found to be equal to 10 hours. 

The relaxation time for transition with frequency cjd ; which involves a simultaneous electron- 
nuclear spin flip-flop, at T = 1.25 K, C ~ 10 16 cm -3 and B ~ 0.3 T was rip ~ 30 hours S> tiib, 7 ][c- 

The extremely long relaxation times of the electron and nuclear spins imply that the required 
initializing of nuclear quantum states (full nuclear nonequilibrium polarizations) can be attained by 
deep cooling of short duration of only nuclear spin system to T\ < 1 mK without deep cooling of the 
lattice. There is the possibility to reach it at the indirect cooling of nuclear spin system by means of 
dynamic nuclear spin polarization techniaues\19\. 

One such method of dynamic nuclear spin polarization for donor atoms is based on the saturation 
by the microwave pumping of the forbidden transition (frequency uj$ in Fig. 3), that is designated as 
the Abragam's solid state effect^ IT9] . 

Let us consider this effect as applied to the ensemble of 31 P atoms. The polarization of electrons 
Ps = 2 (S z ) and of nuclei Pi = 2 (J z ) may be for the sake of simplicity expressed as 

P s = p(l,l)+p(l,0)-p(l,-l)-p(0,0), 

Pi = p(l,l)+p(0,0)-p(l,0)-p(l,-l), (28) 

where p(F, m-p) are the populations of states \F, m^) (Fig. 3). They also fulfill the requirement 

p(l,l)+p(l,0)+p(l,-l)+p(0,0) = 1. (29) 

The rate equations for the populations are (it is assumed, that the relaxation rates for transitions at 
frequencies ui^ are equal to P||a) : 

dp(0,0)/dt = (p(l,l)-p(0,0)r B )/r||B + (p(l ! 0)-p(0 J 0)r D )/r| |D + 

+(p(l,-l) -p(0,0)r+)/7|| A , 
d P (l,-l)/dt = ( P (l,0)-p(l,-l)r c )/T||c + (p(l,l)-p(l ) -l))-W r e + 

+(p(0,0)r+-p(l,-l))/2]| A , 
dp(l,0)/dt = (p(l ) -l)r c -p(l,0))/T|| C +p(0,0)r D -p(l,0))/T|| D + 

+ (p(l,l)-p(l,0)r A )/2|| A , 

dp(l,l)/dt = (p(0,0)r B -p(l,l))/T|| B + (p(l,-l)-p(l,l))-^e + 

+(p(l,0)r A -p(l,l))/7]| A , (30) 

where parameters rB,c,D,A = exp(— ^cjb,c,d,a/^P) are ratio of rates for an up and down thermal 
transitions. For values hujB,c,T>/kT ^> 1, fiuij^/kT <C 1 (T < 0.1K) there are the thermal electron 
Pg w — 1 and nuclear Pio = huo^/kT <C 1 polarizations. 

Let us assume next that the rate of induced forbidden electron transitions |1,1) => |1,— 1) at 
frequency u>s, that is Ws and electron longitudinal relaxation times satisfy the conditions: 

^ST 1 < T ||B ~ T\\C < ^IIA^HD, 71|S, ( 31 ) 
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where rip, rug are the longitudinal relaxation times of electron spins for forbidden transition. Hereafter 
we shall write 

d P (o,o)/dt = p(i,i)/7 1 | B + (p(i,-i)-p(o,o))/r|| A , 

d P (l, ~l)/dt = p(l, 0)/T|| B + ((p(l, 1) - p(l, -1)) • W c + (p(0, 0) - p(l, -1))/3]|A, 

dp(l,0)/dt = -p(l,0)/T|| B + (p(l,l)-p(l,0))/I|| A , 

dp(l,l)M = -p(l,l)/T||B + (p(l,-l)-p(l,l))-We + (p(l,0)-p(l,l))/I]| A (32) 
With equations (|28j ) .l|29 |) .(|31j ) we can obtain the rate equations for Pg and -Pi: 

cLPsM = -(Ps + Pi)-W e -(P s + l)/T l]B , 

dPi/dt = -{P s + Pi)-W e -Pi/T u . (33) 

The steady-state saturation condition (W c 3> l/T!i A ) of the transition |1, 1) |1, —1) gives rise to the 
equalization of the populations p(l, 1) = p(l, —1) and to the full nuclear spin polarization respectively 

Pl = -Pg=p(0,0) = l. (34) 

It is obvious that this state is equivalent to the state with nuclear spin temperature T\ < huop^jk ~ 
1(T 3 K. 

Let us estimate finally the needed microwave power for saturation. The rate of external microwave 
field that induces forbidden electron transitions Ws differs from the rate of allowed transitions with a 
flip of only one nuclear spin W by small factor that is proportional to (Bs/B) 2 jS], where for isotrope 
hyperfine interaction B$ is the field due to the dipole-dipole interaction between nuclear and electron 
spins of 31 P atom. Let us write 

W s ~ (7e&mw) 2 • rig/2, (35) 



where (Aajg) 1 / 2 ~ ^r^g 1 is the nonhomogeneous broadened resonance line width for the saturated 
electron transition, t£ s , is effective transverse relaxation time of electron spins, 6 mw is the amplitude 
of microwave field. 

As a result the saturation condition takes the form 

Ws>1/t|| S ,T|| A or W>l/r m , ( lc b mw ) 2 ■ r* ±s r m > 1, (36) 

that is the same form as for allowed transition. By using expression for quality factor Q c of microwave 
cavity 

Qc « ^S&mw • V I /(2 P L P), (37) 

where V r is volume of microwave resonator, P is dissipated power we will obtain the following saturation 
condition 

W e »l/T|| S or ( 7c 6 mw ) 2 Tt s T,| S » 1. (38) 
The dissipated power in cavity for (7 e 6 m w) 2 ■ T ls T ||B = 1 is determined by 



P > u; s V r W(Aui) 1 / 2 /(2 l i Q cl 2 ). (39) 



For example, taking W ~ 1/t|| B ~ 103s- 1 , uj s ~ lOOradGHz, V v ~ 1cm 3 , Q c ~ 1000 and (Acjg) 1 / 2 ~ 
10 s s — 1 [20j as a rough estimate we obtain P > lmW. Notice that this power is applied only during 
the saturation process over the time > W' 1 ~ 1 ms in the act of qubit state initialization. 
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Hence, the initialization of nuclear states may be obtained by using ENDOR technique at the 
lattice temperature of the order of 0.1 K and by this means that the second difficulty of Kane's scheme 
can be overcome. 

Notice here that there is also another possibility of ensemble NMR implementation, which does 
not have the gate system. The selectivity of nuclear resonance frequencies for individual qubit in the 
ensemble of Kane's chains can be achieved, rather than using the A— gate voltage, with the applying 
of the external magnetic field gradients along axis x. For neighboring qubits separated by ~ 20 nm it 
is required dB z /dx ~ lT/cm (that is feasible now), which produce a resonance frequency difference 



There was proposed previously an all silicon NMR quantum computer where qubits are Si nuclear 
spins arranged as chains in a 28 Si matrix and with natural crystals of calcium hydroxyapatite, 
involving one-dimensional hydrogen chains in both cases nuclear resonance frequencies are sep- 
arated by a magnetic field gradient. However, in this case too large field gradient, of the order of 
1 T/jUm is required. 

6 The nuclear spin states decoherence due to hyperfine interaction 
of nuclear and electron spins 

The relaxation of nonequilibrium state of the nuclear spin system represented by the product of 
independent (nonentangled) one-qubit states, owing to the interaction with isotropic environment, 
shows two processes. One is a slow establishment of equilibrium state associated with dissipation of 
energy. For it the diagonal elements of density matrix decay with characteristic longitudinal (spin- 
lattice) relaxation time Tji . The decay of non-diagonal matrix elements called decoherence of quantum 
states is characterized by a decoherence time or transverse (spin-spin) relaxation time T±. The 
longitudinal relaxation times Tji in the case of nuclear spin of 31 P atoms as qubits is defined mainly 
by thermal modulation of qubit resonance frequency accompanied by spin flips. It is usual that for 
solids Tj. <CI]|. 

The internal adiabatic decoherence mechanisms due to a random modulation of qubit resonance 
frequency, produced by local fluctuating magnetic fields without spin flips. These fields are determined 
by secular parts of interactions of nuclear spins with electron spin of the basic phosphorus atoms, with 
impurity paramagnetic atoms and also with nuclear spins of impurity atoms. We have named this 
mechanism as internal. It seem to be the leading one. 

The modulation of nuclear spin resonance frequency Acv(t), which is determined by the secular 
part of hyperfine interaction, may be written as 



where A(t) = Ao+AA(t), AA(t) is the modulation of hyperfine interaction constant, Aq = 725 radMHz. 
The influence of gate voltage noise on this frequency modulation was studied in 1101 I23| and it is 
not treated here (external decoherence process). 

Another (internal) modulation mechanism of A(t) is the interaction of donor atoms with acoustic 
phonons. It is our belief that for very low temperature this mechanism is not essential |24j . 

Let us consider now the first term in ()40|) . We shall follow the semiclassical model of adiabatic 
decoherence of one-qubit state (Appendix A. 2). The correlation function of frequency modulation 
Aujo,(t) = Ao(S z (t) — (S z )) is determined by the fluctuations of electron spin polarization and depends 
on electron resonance frequency us, longitudinal t\ (hours) and transverse T2 relaxation times. In 
adiabatic case ws = 7s -B > 1/t2 3> 1/ri and we will obtain: 



100 Hz. 



Au(t) = A(t)S z (t) - A (S z ) « A (S z (t) - <5 Z » - AA(t) (S z ) 



(40) 




(41) 



where 



( 



Aul) = Al ((S 2 Z ) - (S z ) 2 ) = A 2 (l - tanh 2 ( 7s ^/2A ; T))/4. 



(42) 
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Now, according to (|74jl. we obtain 

r(t) = (Awl) if (i/n - 1 + expH/ri)). (43) 

For n ps 10 4 s and t ~ T d = 1 s, 1< (Aw§) r x 2 < (ri/T d ) 2 we have the non-Markovian random process 
(slow dampening fluctuations). In this case 

T(t) = (Au 2 s )t 2 /2 (44) 

and the effective decoherence time can by estimated from T d ~ (Acjg) 1//2 . 

The necessary value of decoherence time for the NMR quantum computer clock time ~ 10 -4 s 
should not exceed several seconds. Therefore, let us write the requirement for ^ c hB/kT S> 1 in the 
form 

1/T d 2 ~ A 2 Q {\ - tanh. 2 (-/ e hB/2kT))/A « 2A% exp(--f e hB /kT) < Is" 2 , (45) 

from which we find that the decoherence suppression will be achieved only at sufficiently large B/T > 
30T/K. It corresponds to B = 2T for lattice temperatures T < 0.06 K. 



7 The adiabatic nuclear spin states decoherence due to interaction 
with nuclear spins of impurity atoms. 

The paramagnetic impurity atoms having magnetic moments play also a role of environment for 
nuclear spins in solid state. However decoherence mechanism due to dipole-dipole interaction of their 
magnetic moments with nuclear spins-qubits is suppressed to a large extent at B/T > 30T/K thanks 
to near- full electron spin polarization |24| . 

Another mechanism of one qubit state decoherence is dipole-dipole interaction with not fully 
polarized nuclear spins 1^0 of impurity diamagnetic atoms having concentration Ci^mp. Isotope 29 Si 
with 7i i imp = — 53radMHz/T is one of such atoms. The random fluctuating local field, produced by 
nuclear spins of impurity atoms has the form 

N 

AB a (t) = - J^MnXiAimpfot) - (^,im P (ri))) (46) 

i,/3 



where 



n / x MO TiTi.imp ( x 3r iQ r i/3 \ 

D a A r = <W " ] > ( 47 ) 



i \ i 



T[ is the distance-vector to i-th impurity nuclear spin. 

In this case correlation function of frequency modulation 

(Aws(t)Aufc(0)) = 7i 2 (B z (t)B z (0)) = C 1)imp f J^^A^Mp^^Mp^) ~ (If3M P (r)) 2 )dr (48) 
takes the form 

(Asw(t)Aafc(0)) = (Ao; 2 ) exp(-t/I|| )imp ), (49) 

where 7]|,i m p ~ 10 s IS impurity nuclear spin longitudinal relaxation time of isotope 29 P at low temper- 
ature |THJ . Taking 21i i mp to be much more than T d ~ 1 s, for the determination of allowable impurity 
concentration we obtain equation 

1/T 2 « C Ijimp • { ^ ll ^j n)2 ■ (l " tanh 2 (| 7l)imp | KB/WTj)) , (50) 



12 



where a is minimal distance to impurity nuclear spin which for Si is of the order of 5 • 10 cm . 

For B/T > 30T/K and for spin temperature T\ at which there is near-full polarization of nuclear 
spins 

iji^phB/mi > i (51) 

or for T\ < 0.8 mK, we will obtain that the allowed concentration of the isotope 29 Si is 

Cl,im P % < 4.5 • 10~ 2 %. (52) 

This value can be increased due to the further decrease of nuclear spin temperature T\. For comparison, 
natural abundance of isotope 29 Si in natural silicon is 4.7%. At present the realized degree of cleaning 
28 Si is 99.98%, which does not fully suit for our purposes yet. 



8 Adiabatic decoherence of entangled two qubit states 



In the processes of input of information and logic operation performance some nonentangled initializing 
states of quantum register become entangled. The adiabatic process of transverse relaxation may be 
also the main decoherence mechanism of coherent entangled quantum states. 

As a simple example let us consider here the adiabatic decoherence of the pure fully entangled two 
qubit triplet state EPR-type |V>epr) = \/l/2 (| j [) + \[ j)) with the zeroth projection of the total 
spin on z-axis, which has density matrix 



Pepr = \iPepr) (V'eprI 



( \ 

110 

110 

\ J 



(53) 



The action of the environment on qubit states will be described quasiclassically as correlated random 
modulation of the qubits resonance frequencies A<jj\^(t) and of indirect spin-spin interaction parameter 
Au(t) = AI\{t)/2. The secular part of Hamiltonian for interaction with the environment is represented 
by 



H(i) = -Awi(t)(cri z ®l)/2- Acj 2 (t)(l®cT2 Z )/2 + 
+Au I (t)(a lz ® <T 2 z)/2, 



(54) 



where 0"i Z) 2 z are Pauli matrix. 

The density matrix Q53|) under the action of random field in rotating frame with resonance frequency 
loq is described by expression 



PEPR(i) = V(ty l p EPR V(t). 

In the considered case unitary matrix 4x4 U(t) is ( = Jo Aui^dt, <pi(t) = Jq Auj\dt 

U(t) = (cos((^i(t)/2)l+isin((^i(t)/2)cJi z ) ® (cos(^ 2 (t)/2)l+i sin((^ 2 (t)/2)cT2 z ■ 
•(cos<£i(£)(l ® 1) + isunpi{t){ai z <g> cr 2z )). 



(55) 



(56) 



For perturbed density matrix we obtain 







PeprW = ^ 



( 

1 eM-i(Mt) ~ Mt))) 

exp(i(y»i(t)-^2(*))) 1 
\ 






o J 



(57) 



We see that the modulation of spin-spin interaction has no effect on density matrix of triplet EPR- 
state. 
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Let us assume now that the random phases <pi >2 (t) have mean value (<pi %(t)) = and belong to 
the reduced statistical ensemble, which is described by two-dimension Gaussian distribution: 



1 



2ira l (t)a 2 (t)^{l 
f 1 



pUt)) 



■ exp 



<pl(t) 2p 12 (i)p 1 (% 2 (i) <p 2 2 (t) 



2(1 -p( 2 (t)) \a((t) 



+ 



(58) 



Here 



are the variances and 



< 2 (*) = (^, 2 (<)) = 2 [ (t - T)h, 2 (r)dT 



(<Pi(t)Mt)) _ 2j^t-r)f 12 (r)dr 



o-i(t)a 2 (t) 



ox{t)a 2 (t) 



where 



/ 1)2 (r) = (Awi i2 (r)Au;i,2(0)) , /i 2 (r) = (Aa; 1 (r)Ao; 2 (0)) 
The normalized mutual correlation function p\ 2 (t) takes values in interval 

< pi 2 (t) < 1. 
After an averaging (|5T)) with (|o^|) we have 



(59) 



(60) 



<Pepr(*)) 



/ 
1 

exp(-r(t)) 

V o o 



\ 
exp(-r(£)) 

1 
J 



(61) 



where 



exp(-r(t)) 



dtpi 



dip 2 w(<pi,<p 2 ) exp(±i(<pi - ip 2 )) 



exp{-(al(t) - 2a 1 (t)a 2 (t)p 12 (t) + <r|(t)}. 



(62) 



In the absence of random field correlation pi 2 (t) = decrement T(t) is equal to the sum of decrements 
of two one qubit states: 



r(t) = (^(i) + ( 7 2 2 (i))/2 = 2r 1 (t). 



(63) 



In case of maximum correlation p\ 2 (t) = 1 and (fi(t) = (p 2 (t) (the same mode acts on both qubits) 
adiabatic decoherence disappears. Analogous properties have the singlet EPR state. 

We see here, that decoherence of interacted qubits states may differ essentially from one qubit 
decoherence. Under the action of fully correlated random fields the coherence of two mentioned 
entangled states is not violate and they may be considered as the basis of decoherence-free substrate 
for logical qubits coding. Clearly the pure nonfully entangled states \ip) = — a \ j J.) + ^fa \ it)) 
have no such properties. 

Adiabatic decoherence of other two qubits fully entangled quantum Bell states \ip) = \f\j2 (| f j) ± 
III)) under the action of fully correlated random fields with pi 2 (t) = 1 now does not disappear. Its 
decrement is now equal to T(t) = (<7i(i) + a 2 (t)) 2 /2 = 2a\ = 4Ti(t), that is four times larger than 
for one qubit decoherence. 
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9 An antiferromagnet-based ensemble NMR quantum computer of 
cellular-automaton type 



For implementation of ensemble silicon quantum computer, operating on cellular-automaton principle, 
it may be usable the previously considered ensemble of long chains of donor atoms 31 P disposed in 
silicon, but free of the A and J gates. 

If exchange interaction constant (it is here positive) for localized electronic spins of 31 P along the 
chain is more than Zeeman energy J(l) 3> ^ehBr ~ 6.5 • 10 -23 J, that corresponds to the distance 
between donors l x ~ 20 nm JO]; that is the electron critical temperature (Neel temperature) will be 
Tns ~ 4K. and the lattice temperature is well below the critical temperature for electron ordering 
T < Tns ~ J{l)/k (k = 1.38 • 10~ 23 J/K is the Boltzmann constant), than the one-dimensional 
antiferromagnetically ordered ground state of electronic spins can be produced. 

Due to hyperfme interaction nuclear spins will be oriented according to the electronic spin direction 
in the resultant field and can form array with the alternating orientation of nuclear spins. At magnetic 
fields B < A/^/ih ~ 3.5 T and at spin temperatures T ~ 10~ 3 K the nuclear spins 31 P will form a 
periodic ground state array of ABAB. . .type: Till • • •> where j marks the ground state of nuclear 
spin in an A-site and [ is the ground state of nuclear spin in a B-site with almost 100% opposite 
orientation {tox,B/kT > 1). That is the distinct nuclear spins will be in initialized ground state. 

Notice that the using of dynamic methods makes possible the high orientation of nuclear spins also 
at larger lattice temperatures, this state will be the long-lived nonequilibrium nuclear spin state. 

The nuclear resonant frequencies wa,b of neighboring nuclear spins are different for each of the 
magnetic one-dimensional subarray A and B in the chain as they depend on the states of neighboring 
spins. We will take it here in the simple form |14l 1251 l2l)] : 

wa,b = oj(m < + m>) » \jih~B ± A/2 - I n ■ (m< + m>)| /%, (64) 

where I n is constant of two neighboring nuclear indirect spin-spin interaction, m < = ±1/2 and m> 
= ±1/2 are the magnetic quantum numbers for the left and right nuclear spins. The nonsecular part 
of nuclear- nuclear interaction is neglected here taking into account that jihB, A/2 S> I n - The differ- 
ence between nuclear resonant frequencies for the distinct neighboring spin orientations is Aoji/2tt ~ 
I n /27rh ~ 0.5 MHz, whereas the resonant nuclear frequencies are cja,b/27t ~ A/Airfi ~ 120 MHz. 

For the organization of logic operations let us use the addressing to spin states, similar to the 
scheme put forward in [22]. Each nuclear spin in A-site of this scheme has two internal eigenstates - 
ground | f ) and excited |J|) and in B-site -\\) and |ff) accordingly. 

We take into account that the life time of excited states (the longitudinal nuclear spin relaxation 
time Tji ) at low temperatures is very long. Each logic qubit of quantum information in this state will 
be encoded here, similar to [2*7]. by the states of four physical spin-qubits: the logical qubit basis state 
"0" will be encoded by unit l-U-ffT l)> while the state "1" will be encoded by | j j -IJ-ff ). It is important 
here that the resonance frequencies of nuclear spins depend on neighboring spin states. Both logical 
states have two excited spin states and zero projection of total nuclear spin. 

Notice that a random inversion of only one spin will result in degradation of the qubit state. But 
to form the rough error, for example, of "0" "1" type in the coding of stored quantum information 
it is essential to invert four spins simultaneously. Therefore, it may be concluded that the considered 
way of qubit coding ensures a better fault-tolerance with respect to this type of errors. 

The input and output of the information in the array of ground states spins could be performed at 
the ends of the array, where the nuclear spins (say in A-site at the left end) have only one neighboring 
spin and resonant frequency wa(- 1/2) (m< + m> = —1/2). The corresponding selective resonance 
7r A,-i/2 — pulse inverts only one nuclear spin (in A-site) at the end of array and doesn't influence 
other ones. Then the new selective 7Tb,o — pulse will invert next nuclear spin (in B-site), which has 
the opposite orientation of ground and excited neighbor nuclear spin (m< + m> = in A-site) and 
consequently the new resonant frequency, distinguished from the frequency of spins with the neighbor 
nuclear spin in the same ground states (m< + m> = 1) (See Table). 
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Table. The tt— pulses for spins in A- and B-sites 



Neighbor spin states A— site 


A I 


I A I 


I A^ 


it A I 


H At 


Resonance frequency 


i/a(-1/2) 


^a(-I) 


^a(O) 


^a(O) 




7T— pulses 


7TA-1/2 


7TA.-1 


TTA.O 




7TA.1 


Neighbor spin states B— site 


B I 


t B t 


T B JJ- 


4 B t 


4 B 4 


Resonance frequency 






^b(O) 


^b(O) 


fB(-l) 


7T— pulses 


7TB.1/2 


7TB, 1 


7TB,0 


7I"B,0 


7TB, -1 



Thus the logical qubit state " 0" , that is | JJ- f J, ) , is formed in the following way (the pulses act on 
underlined spins from the left): 

"0" 

TITI •••VT A ,-l/2 4ITIT •••VTB.O MTIt ••• 

ABAB AbAB abab 

As the ports (noted in Fig. 5 by cross 'x') for input and output of the information in the array of 
ground states spin-qubits can be also dopant nuclei D at the certain place of the array with distinct 
resonant frequency, defects or local gates that modify the resonant frequency of the nearest nuclear 
spin in the array. 

Electron spins S 

Nuclear spins / J I x I 

port 

Fig. 5. Scheme of electron and nuclear spins ordering. 

Starting from the perfectly initialized states inputting the information can be performed by setting 
the dopant D-spin to a desired state by means of RF-pulse at its resonant frequency. The nuclear spin 
state of spin nearest to the dopant spin is created by SWAP operation. After the required information 
is loaded, D-spin is reset to the ground state. Upon completion of computation, the state of any spin 
can be measured by moving it to the A-site nearest to D, then swapping A D and finally measuring 
the state of D-spin. 

For the implementation of quantum operations on logic qubits we will also introduce, one auxiliary 

control unit (CU), which is represented here by six physical spin states in the pattern ft-lj- 1 TftJJ-- The 
CU exists only in one place along the array and is separated from logical qubits by odd number 
of spacer spins. The applying of corresponding SWAP sequence of pulses CU leads to putting the 
interaction of CU with one and two logical qubits and performing on them one- and two-qubit quantum 
operations (for more details, see [T H 125 } T26 ) . 

In the case of a large enough ensemble of in parallel acting chains the states may be measured by 
NMR methods. For the increasing of logic qubit number in single 'molecule' to > 10 3 , that falls on 
one port, two- and three-dimensional structures with antiferromagnetic chess-type ordering of electron 
spins may be used. The corresponding ordering will also be for nuclear spins (Fig. 6). 
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port _>| / x |<_ 

J t | i | t | t | t | — 
t i t i t J t J t J t -i- 
t i t J t J t J t J t 

Fig. 6. Scheme of two dimensional chess-type ordering of initialized nuclear spins. 

If, for example, the number of spins-qubits, falls on one port in linear chain is, say, L ~ 30, so in 
two-dimension case their number will be L = 900. Ensemble that is composed of N ~ 10 5 in parallel 
acting such plane 'artificial molecules' permit to provide input and output of information through the 
standard NMR techniques. The using of more sensitive ENDOR techniques has particular meaning 
when it may be combined with the techniques of dynamic polarization (solid state effect). 

Structures with two and three-dimensional antiferromagnetic order may be found perhaps among 
the natural rare earth or transition element dielectric compounds. 

There are the rare earth compounds of thulium stable isotope 169 Tm, that has nuclear spin / = 1/2, 
#N = 0.458 and makes up 100% of abundance, with stable spinless isotopes of other elements Tm. 
They can be possible: Trr^Os, TmSi2, TmGe2, TmSe. The natural elements O, Si, Ge and Se have, 
accordingly, nuclear spin containing isotopes (in brackets the isotope abundance is shown) 17 1 = 5/2 
(0.04%), 29 Si I = 1/2 (4.7%), 73 Ge I = 9/2 (7.76%), 77 Se / = 1/2 (7.78%). The choosing the needed 
compounds requires further detailed theoretical and experimental investigations. 

It may be also considered, as variant for ensemble NMR quantum computer, organic dielectric crys- 
tals, containing quasi-one-dimensional chains, such as antiferromagnetically ordered chains of poly- 
acetylene with only proton nuclei, for qubits: 

H H H 

c c c 



12 c c c 



H H 



The advantages of ensemble quantum cellular automaton in comparison with the above considered 
ensemble variant with strip gates are as follows: 

a) The system of the control gates is absent, what essentially simplifies the production of computer 
structure and eliminates one of the important sources of decoherence. 

b) The coding of logic qubits into four physical qubits gives a higher degree of fault-tolerance in 
logic operations. 

It follows that the going to ensemble quantum cellular automaton permits to overcome the third and 
fourth difficulties of Kane's scheme. 

The chief disadvantage of cellular automaton scheme is the relative complexity of logic operation 
performance. 
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Conclusion 



1. The line of the large-scale ensemble NMR quantum computer development has certain advantage 
over Kane's scheme. It consists in the possibility of employment of the standard NMR technique 
for the measurement of quantum states at output of computer, like in the liquid prototype. 

2. For the initialization of nuclear spin states at temperature T ~ O.IK methods of dynamic 
polarization may be proposed. 

3. Analysis of proposed planar structure of ensemble silicon computer shows the possibility of 
realization of large-scale NMR quantum computer for ensemble component number N ~ 10 5 . 

4. The main reasons for the internal decoherence of one qubit states are the modulation of resonance 
qubit frequency due to hyperfine interaction with fluctuating electron spin and due to interaction 
with randomly distributed impurity diamagnetic atoms containing nuclear spins. 

5. Analysis of different feasible ways for obtaining decoherence times large enough shows that the 
values, needed to perform the required for large-scaled computations number of quantum logic 
operations ~ 10 5 , can be achieved. 

6. The implementation of cellular automaton principle permits to abandon the realization regular 
nanostructure in the form of gate chains. 

Appendix 

Al. Signal NMR for discrete ensemble of nuclear spins 

Let us consider here the sample that involves N = nL ■ pNo nuclear spin-qubits arranged in the plane 
z = of the silicon plate at regular intervals along the strips (Fig. 4) . The spins in chains under strip 
at resonance in each block are oriented along x axis (solenoid axis) and separated at intervals of L. 
The read out NMR signal is 



K f x / 2 f 
\V max \ = Qo;a— / / B xmax (x,y,z)dydz 
A J -xii J A 



k r x / 2 

-X/2 

where X = nL » L is roughly the length of solenoid and 



dx, (65) 



- ( ' 1 " «* n h n l<* - W + (v - W + = 2 P 2 (w) 

is the peak magnetic field produced by resonant spins in solenoid. For simplicity it is suggested that 
n,p, No, L S> 1 are the even numbers. We assume that the area of coil turns is A = D ■ 5 (D = l y -pNo, 
<5«D). 

For summation over n\ and p\ we have used the Poisson summation formula, namely, 



pNo/2 



nNa 00 r D / 2 

E /(W = E 7rE/ /(6exp(^/Z y R, (67) 



Pi =-pN /2 v=- x 

by omitting the oscillated terms with v / 0: 



jlh 1 f x / 2 n f x l 2 

l^maxl = UoQKuJA — — / dx— / df]- 

8tt X J -X/2 X J-X/2 

[ 5/2 d7 P^o f D / 2 f D ' 2 ~2(x - V ) 2 + (y- Q 2 + z 2 
Jo D Lb 12 J-D/2 [(x - n) 2 + (y - £) 2 + z 2 }*/ 2 ^ 



(68) 



18 



Taking into account D 2 S> z 2 , upon integrating (|67j) over y, £, we obtain 
|T/ max | = ^^——-J^dx-J^dr,- 

5/2 dz ( ^ ~ ~ ^ + ^ ^ (69) 

\[(x - rf) 2 + z 2 ] 2 [D 2 + (x - rf) 2 ] 1 / 2 [(x - V ) 2 + z 2 } 3 / 2 J ' 1 ; 

By integrating now over z and preserving only the major logarithmically increasing for 2\x — r/\ jb => 
term, we obtain 

-fihnpNo 1 f x / 2 , 1 [ x ' 2 



Vmax ~ hoQKAuja- — — — / cte— / dl] log 



X/2 X J-x/2 \x - rj\ 

N X X 

= (fi /4)QKAu; A - 11 h-—log TT , (70) 

where V s = AX, N = npN , e = 2.718. . . . 

We see that expression (A. 7) is distinguished from (j3J) by non-essential factor ^log^j, that is 
of the order of several ones. 



A2. Semiclassical model of adiabatic decoherence of one-qubit state 



We will consider a long-lived non-equilibrium qubit state when diagonal elements of density matrix 
may be treated as a constant. 

The random modulation of resonance frequency Au;(t) that causes the dephasing of a qubit state 
is determined by the random phase shifts 



<P(t) 



iu(t)dt. 



(71) 



The one-qubit density matrix of pure state in rotating frame with non perturbed resonance circular 
frequency will be 



pit) = 1/2 



1 + P z P_ exp(zp(t)) 

P+exp(-vp(t)) \-P m 



(72) 



P2 + P 2 +P 2 



where P± = P x ± iP y , P x , P y ,P z are Bloch vector components of length P - 

By treating the resonance frequency modulation as Gaussian random process after averaging 1)72(1 
over phase distribution with {f(t)) = we obtain 



<p(t)> = 1/2 



1 + P z P_ exp(-r(t)) 

P+exp(-r(t)) 1-P Z 



where 



r(t) = 1/2 



Au(t)dt 



(t - t) (Aw(t)Aw(O)) dr, 



(73) 



(74) 



f(t) = (Auj(t)Auj(0)) is the frequency correlation function of a random process, which is characterized 
by variance (Au;(0) 2 ) and correlation time tq such that for t > tq (Auj(t)Auj(0)) 0. For T(t) > 
the averaged density matrix presents a mixed quantum state with two non-zero eigen states 



1/2 ■ 1 ± Jl - (P 2 + P 2 )(l - exp(-2r(t)) 



(75) 

and the populations of states p± = 1/2(1 ± P z (0)) at T(t) oo. 

Thus, the adiabatic decoherence problem is reduced to the determination of the function T(t) or 
the correlation function of random frequency modulation. 

In the case of an ensemble quantum register there is a need to average the one-qubit density matrix 
and correlation function over ensemble of independent equivalent spins-qubits. 
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